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We discuss the possibility of spin glass order in the vicinity of the unexpected metallic state of the 
two-dimensional electron gas in zero applied magnetic field. An average ferromagnetic moment may 
also be present, and the spin glass order then resides in the plane orthogonal to the ferromagnetic 
moment. We argue that a quantum transition involving the destruction of the spin glass order in 
an applied in-plane magnetic field offers a natural explanation of some features of recent magneto- 
conductance measurements. We present a quantum field theory for such a transition and compute 
its mean field properties. 



The pioneering experiments of Kravchenko, Pudalov, 
Sarachik, and others Q have demonstrated that the low 
density electron gas in two dimensions (the 2DEG), in 
the absence of an applied magnetic field, exhibits some 
fascinating strong correlation physics. The initial indi- 
cations of this new regime came from an unexpected de- 
crease in the resistivity as the temperature was lowered 
on the electrons in a silicon inversion layer. Experimen- 
tally it is of course impossible to definitively rule out that 
there will be an eventual upturn of this resistivity at some 
exponentially small temperature, and the that the true 
ground state is always an insulator. Nevertheless, the 
single-electron-transistor measurements of Hani et al. ||] 
show that there is a remarkable qualitative change in the 
fluctuations of the local compressibility at the same den- 
sities, and there must be a corresponding transition in 
the quantum ground state. A theoretical description of 
this correlation-induced transformation remains an im- 
portant theoretical challenge. 

Useful guidance is obtained from numerical studies of 
the disorder-free 2DEG at the same densities. These 
indicate that the system is not too far from a Wigner 
crystallization transition, and that non-trivial spin cor- 
relations are present ||. A study of an effective spin 
Hamiltonian |||5| for the Wigner crystal ground state 
shows instabilities to ferromagnetism || and to a spin 
singlet ground state which is an attractive candidate for 
fractionalization. 

Turning to the disordered 2DEG, we can reasonably 
expect that there will be regions of the sample where the 
system is locally ferromagnetic However, in other re- 
gions the ring exchange terms may prefer a spin singlet 
state, and these should couple the different ferromagnetic 
regions so that they align in different directions. On these 
grounds, we propose here that the ground state has spin 
glass order in the vicinity of the transition in the con- 
ductivity, possibly with a concomitant average ferromag- 
netic moment; if the ferromagnetic moment is non-zero, 
the spin glass order resides in a plane (in spin space) or- 
thogonal to the direction of the moment. We expect that 
the spin glass character is stronger on the metallic side, 
while the ferromagnetic character of the ground state be- 



comes stronger as the density of electrons is reduced, and 
this is partially responsible for the decrease in the con- 
ductivity. Other arguments for the proximity of the spin 
glass order, which were based directly on the theory of 
the disordered metal, were presented in ||. 

Our proposal is also motivated by recent experimental 
studies of the magnetization of the 2DEG |p| -H . In par- 
ticular, the measurements of Vitkalov et al. |9| in an in- 
plane ("parallel") magnetic field, H, indicate that there 
is a well-defined critical field H a , dependent on the den- 
sity of the 2DEG, at which the magnetoconductance dis- 
plays a kink-like feature, and which they identify with a 
quantum phase transition. For a clean 2DEG it is natu- 
ral to interpret H a as a saturation field, above which all 
electrons have their spin polarized in the direction of the 
applied field. In a disordered 2DEG there is no reason for 
all regions of the sample to respond in the same manner, 
and there will always be local configurations in which the 
electrons prefer to form low spin clusters at these rela- 
tively weak fields. Naively, one would then conclude that 
the magnetization of a disordered 2DEG should evolve 
smoothly as a function of the applied field, without non- 
analytic structure at any critical field H a . A simple way 
out of this dilemma is to search for some other order pa- 
rameter associated with the preserved symmetry of spin 
rotations about the direction of the applied field. Spin 
glass order in the plane perpendicular to the magnetic 
field is our proposed candidate: this order is present for 
H < H a , but vanishes for H > H a when the polar- 
ization of many spins along the direction of the applied 
field reduces their "canting" enough to disrupt temporal 
memory of their orientation in the orthogonal plane. 

(The discussion of this present paper will restrict at- 
tention to the case where H a > 0. Experimentally, it 
is known that H a decreases as the density of electrons 
decreases, and a possible reason for this is the increase 
in the spontaneous ferromagnetic moment which acts as 
an effective magnetic field. It is possible that there is a 
density at which H a reaches zero: this corresponds to a 
transition at H = 0, driven by tuning the density, be- 
tween a ferromagnetic state with spin glass order in the 
transverse plane and an ordinary ferromagnet. The the- 
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ory for this density-tuned transition will be similar to 
the field-tuned transition described below, but will not 
be presented here: it follows by combining the results 
below with those of 

The presence of a critical point at H = H a also leads to 
a natural explanation of another puzzling feature of the 
data. For the density at which H a was small, it was found 
that the characteristic iJ-width of the magnetoconduc- 
tance feature scaled roughly as the absolute temperature, 
~ UbT (we have absorbed a factor of gps in the defini- 
tion of H, where \xb is the Bohr magneton and g is the 
g-factor of electrons in the silicon conduction band.) A 
trivial system whose thermodynamics depends only on 
H/T are isolated free spins. However given the large ex- 
change interaction energies between the electrons, it ap- 
pears extremely unlikely that there is a sufficient density 
of isolated free moments to lead to a significant change 
in the magnetoconductance. Indeed there is a 'catch- 
22' here: if the moments are really isolated enough to 
behave as free spins, their coupling to the itinerant elec- 
trons is weak and they have a negligible effect on the 
conductance. The presence of a spin-glass quantum crit- 
ical point at H = H a offers an alternative route to ob- 
taining a characteristic field of order T. Indeed it can 
be argued on rather general grounds that if the quantum 
critical point obeys the hyperscaling property, and if H 
couples to a conserved total spin, then the characteristic 
field scale will be of order ksT [ p"5| . 

It is appropriate to now mention very interesting recent 
observations of glassy behavior, persisting in the metallic 
phase, by Bogdanovich and Popovic [ fl6| |, which appear 
to be consistent with our proposal. These authors made 
the connection to glassy behavior in models of charge 
transport of spinless electrons jl7j], while here we will 
focus only on the spin degrees of freedom. 

We will now propose a quantum field theory for the 
transition at H = H a . We will follow general approach 
to quantum spin glasses developed in ]12j and reviewed 
in Chapter 16 of |lq |. In the presence of an applied mag- 
netic field, there is only a U(l) symmetry of spin rotations 
about the direction of the applied field, and we are inter- 
ested in singularity in the dependence of the conserved 
"charge" of this symmetry (the magnetization) on the 
applied field. A general theory for such transitions was 
presented in and in Chapter 11 of [Q. Here we will 
extend this theory to random systems with spin glass 
order by combining it with the methods of |]l2f . 

Let us assume that the field H is applied in the z di- 
rection, and let S a (r, r) with a — x, y be the component 
of the electron spin density in the orthogonal x,y plane at 
the spatial point r and at imaginary time r. We use the 
standard replica method to treat the quenched disorder: 
consequently, we introduce replica indices a = 1, . . . ,n, 
and the spin density S aa (r, r). The quantum theory for 
spin glass order is expressed in terms of the order param- 



eter functional 7*1, t 2 ) which is 



Q%( r >n,T2) ~ S aa (r,T 1 )S l3b (r,T 2 ). 



(1) 



By applying the methods developed in |12].|18|] to spin 
glass order in the x,y plane in the presence of an applied 
magnetic field we obtain the following low order terms in 
the effective action for Q^g(r, n, t 2 ): 
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where d = 2 is the spatial dimensionality, summation of 
the spin indices a, j3, p over x, y is implied, e a p is the 
antisymmetric tensor, is the mean-field value of the 
critical field (which will be renormalized by fluctuations), 
k, u, v are coupling constants, and we have omitted an 
additional quadratic term which has no influence on the 
mean-field theory, but does play an important role in 
the violation of hyperscaling in the perturbative fluctua- 
tions |L2| . In a metal with an appreciable contribution of 
low energy spin excitations associated with the particle- 
hole continuum, there will also be additional dissipative 
term in the action, as discussed in [ p0[ : in the present 
situtation with strong local correlations we believe this 
is unlikely to be the case, and so have omitted such a 
term above. 

An important property of (||) is that the field H cou- 
ples to a conserved U(l) charge: consequently changes 
in H can be absorbed by a time-dependent gauge trans- 
formation which is equivalent to a transformation into a 
"rotating reference frame" [|l5|jl8| . More specifically, if 
we generalize the action to a time-dependent field H(t), 
then it is invariant under the infinitesimal transformation 

Q* a 0(r,n,T 2 ) -> Q a a a l3 {r,T 1 ,T 2 ) - e a7 <f>(n)Q 7 p(r,n,T 2 ) 

- ep~ / 4>(T2)Q ai (r,n,T2) 
H{T)^H{r)-id T <i>{T), (3) 

where 4>{r) is infinitesimal. If the critical point at H = 
H a satisfies strong hyperscaling properties, then (|^) im- 
plies that the free energy density T sg obeys 
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where the exponents /i = l+d/z, ip = 1, z is the dynamic 
critical exponent, and $ is a scaling function. Taking the 
H derivative of M) we obtain the magnetization 



M = Mn - T^-v®' 



H - H a 



(5) 



where Mq is the background ferromagnetic magnetiza- 
tion present for H 3> H a . The onset of spin glass order 
causes a decrease in this magnetization as H is lowered. 
A related scaling form should also hold for the magneto- 
conductance. 

We illustrate the above behavior of the magnetization 



by a simple mean field analysis of S sg in 



Unfortu- 



nately these mean field results do not satisfy hyperscaling 
properties appropriate to any value of d, and the reasons 
for this are similar to those discussed at length in 
for other quantum spin glasses. The actual situation for 
realistic spin glasses remains an open problem, as an anal- 
ysis of fluctuations about the mean field solutions leads 
to a runaway flow to strong coupling. If our proposal 
is indeed the correct explanation for the experimental 
data |^] , then the hyperscaling prediction tp = 1 must be 
valid for d = 2. 

At the mean-field saddle-point, we may make the fol- 
lowing ansatz for Q a b p(r, n, T2) based upon the require- 
ments of translational invariance in space and time (we 
henceforth set k E = 1): 



iuj n (ri—T2) 



(G) 



where uj n is a Matsubara frequency, and 



Inserting into (g) we obtain the following expres- 
sion for the free energy density after taking the limit 
n -> 0: 

Ts g = T + - ^{(-iw„ +H- H a )D+_{uj n ) 
+ (iu n + H- H a )D_ + {u n ) + 2q EA (H - 

- Kq EA {D 2 + _(Q) + D 2 _ + (0)} 



Qea 



Qea 
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(9) 



we have replace H® by H a anticipating the mean-field po- 
sition of the critical point. It is now straightforward to 
determine the saddle-point of T sg with respect to vari- 
ations in — (a->„), D — y(u) n ), and q E A- We find two 
classes of solutions describing the paramagnetic and spin 
glasses phases respectively, and we will discuss their prop- 
erties in turn. 

In the paramagnetic phase we have q E A = 0, and 
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D+-(u) n ) = y/—iu)r, 

K 



D- + (LJ n ) = — sj iu> n + A 



(10) 



where the energy A is determined by the solution of 
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H-H a - {u + 2v)Tj2v r - 



(11) 



Dapi^n) = 8 a bD a p(uJ n ) + 5 a p- 



T 



-qEA- 



The first, replica-diagonal term in ^ is the local dy- 
namic spin susceptiblity in the x,y plane, while q E A is 
the Edwards- Anderson spin glass order parameter. For 
simplicity we have assumed a replica-symmetric form for 
the spin glass order — replica symmetry is not broken for 
the terms included in (Q), but is broken when higher 
order terms are included: this is as discussed in Jl2| . 
Notice also that the second term in ([?]) also includes 
a contribution of q E A along the replica diagonal — this 
is unlike the usual procedure in the theory of classical 
spin glasses; the diagonal contribution here accounts for 
the long-time limit of the local spin correlation function. 
Further analysis is simplified by rewriting the Dap in 
"circularly-polarized" components: 



D XX = Dyy = \ (£>+_+£_+) 



D, 



-D 
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D-+). 



(8) 



The condition A > delineates the boundary of the para- 
magnetic phase. At T = 0, the Matsubara summation in 
( |TT| ) becomes a frequency integral which evaluates to zero 
(the expressions in (|l^,|ll|) only aim to capture the sin- 
gular low frequency behavior, and it is assumed that the 
high frequency form is such that contours of frequency 
integration can be freely closed at complex infinity). So 
A = H — H a at T = 0, which demonstrates that the para- 
magnetic phase is stable only for H > H a . For T > 0, 
( JTT| ) can be analyzed using methods that have been dis- 
cussed in some detail in Chapter 15 of |l8| ]; for H — H a 
small, its solution can be written as 



A + (u + 2v)TVA = H -H a + {u + 2v)T 3/2 E 



H_H* 
T 

(12) 



with the function 5 is given by 
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d(y)Vy, (13) 
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where V denotes a principle part and 0(y) is the unit 
step function. The expression (|l^) was derived for 
y > (H > H a ) but has been written in a manner 
which defines it for real y. Despite appearances, the 
function 3(y) is actually analytic at y = 0, and in- 
deed it is analytic at all real values of y; for small y, 
E(y) = C(3/2)/(2VtF) + 0.411958y + . . .. For y < 0, the 
results ([l2{l3) apply in the paramagnetic portion of the 
phase diagram present at T > 0, H < H a (see Fig. |l| 
below). At the critical field, H = H a , ( |l2"| ) predicts 
that A = (u + 2u)S(0)T 3 / 2 ; application of hyperscal- 
ing to this quantum critical region would have implied 
A ~ T x a function of (H — H a )/T, and so it is evi- 
dent that hyperscaling is not obeyed by the mean-field 
theory. Imposing the condition A = in ( |l2| ) deter- 
mines the mean-field boundary of the spin glass phase as 
H = H a — (u+ 2u)S(0)T 3 / 2 at small T: this leads to the 
phase diagram in Fig [j]-the structure of the crossovers 
is very similar to those discussed earlier for other spin 
glasses, and the reader is referred to Chapter 15 of 
for a review. Note that the paramagnetic phase extends 



T 



Quantum 
critical 




FIG. 1. Phase diagram of the quantum spin glass described 
by Ssg in (^). The full line is a phase transition marking the 
mean-field boundary of the spin glass phase-in d = 2 spin 
glass order may only be present at T — 0. The dashed lines 
are crossovers into the quantum-critical region, which survives 
fluctuations in d = 2. With hyperscaling, the characteristic 
energy in the quantum-critical region ~ T; the mean-field 
theory has this energy ~ T 3//2 . 



to H < H a for T > 0, and that the expression (|12|) re- 
mains valid in this regime where y < 0. The free energy 
in the paramagnetic phase is also easily obtained from 
Wrf), and we obtain 



(A-H + H a f _ AT 5 / 2 

K 2 (u + 2l>) 37TK 2 



e (fi+A)/T _ l 
(14) 



The T and H dependence of the magnetization is deter- 
mined by taking a H derivative of ()l4|): we find 



M 



H, 



H 



FIG. 2. Mean field prediction for the magnetization, M, as 
a function of H at T = in (|l|) and @. Upon including 
fluctuations, the singularity at H — H a will survive but the 
magnetization will be H dependent for H > H a . At T > 
the singularity at H = H a will be rounded out on a scale T 
if hyperscaling is obeyed by the critical point. 



up to exponentially small terms, indicating that the mag- 
netization is effectively saturated in the paramagnetic 
phase in this mean-field theory, as indicated in Fig ||. 
Fluctuations will induce a more appreciable variation in 
the magnetization even at T = 0, as strict saturation is 
not possible in a disordered system. 

Finally, we describe the saddle point of (JsJ) in the spin- 
glass phase. Here we find the simple solution 




k(u + 2v) 



(16) 



The last expression identifies the same boundary of the 
spin glass phase (where qEA — 0) as that determined 
above. The free energy can also be computed as before, 
and we obtain 



To 



(H a - H) 2 4T 5 / 2 r(5/2)C(5/2) 
k 2 (u + 2v) 3lTK 2 



(17) 



and the magnetization follows as its H derivative. Now 
we find 



M = Mo 



1(H a - H) 
k 2 (u + 2v) 



as T — ► for H <C H a . (18) 



M = Mo as T -> for H > H a 



(15) 



Comparing with ( |15| ) we see that there is a kink in the 
magnetization at the quantum critical point H = H a , 
as shown in Fig |^. This singularity survives fluctuation 
corrections, even though the mean-field saturation of the 
magnetization in dl5|) does not. 
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This papers has outlined a scenario by which the 2DEG 
can exhibit a quantum phase transition at a critical in- 
plane applied magnetic field H = H a : the transition is 
induced by the destruction of spin-glass order in the plane 
orthogonal to the applied field. We have argued that 
the instabilities of the spin exchange Hamiltonian for the 
ordered Wigner crystal to ferromagnetism and to spin- 
singlet states lend support to the possibility of spin-glass 
order in the disordered 2DEG. It would be interesting the 
compare the H, T, and density dependent data to scaling 
forms like (^,|^)- Theoretical analysis of the transport 
properties of the field-tuned transition at H = H a , and 
also of the density-tuned transition for the case H a = 0, 
will also be of interest. 
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